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Abstract: An interesting alternative to super symmetry (SUSY) for extending physics 
beyond the Standard Model is a model with universal extra dimensions (UED), in which 
the SUSY superpartners are replaced by Kaluza-Klein excitations of the Standard Model 
particles. If new particles are discovered at the LHC, even if their mass spectrum favours 
SUSY or UED, it will be vital to distinguish between their spin assignments in the two 
models as far as possible. We extend the method proposed by Barr [1] to the UED case 
and investigate the angular and charge asymmetries of decay distributions for sample mass 
spectra of both SUSY and UED types. For hierarchical ('SUSY-type') mass spectra there 
is a good chance of distinguishing the spin structures of the two models. However, a mass 
spectrum of the quasi-degenerate type expected in UED would make it difficult to observe 
spin correlations. 
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1. Introduction 

Tiie search for physics beyond the Standard Model (SM) will be a principal objective of 
the main experiments at the CERN Large Hadron Collider (LHC). The leading contender 
amongst theories of new physics is undoubtedly supersymmetry (SUSY), in which the SM 
particles have massive superpartners that differ from them by a half-unit of spin. The 
partners of strongly-interacting partons (the spartons) should be copiously produced once 
their production thresholds are passed. In R-parity conserving SUSY, the spartons typically 
decay into partons and electroweak sparticles, which themselves decay sequentially through 
the emission of SM particles to the stable lightest super symmetric particle (LSP). The LSP 
is usually the lightest neutralino x5, which escapes from the detector unobserved. Thus 
the expected signature of SUSY is the copious production of high-energy jets and leptons 
plus large missing energy. 

If the above signature of SUSY is observed at the LHC, it will be essential to confirm 
as far as possible that the produced objects are indeed superpartners and not some other 
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manifestation of new physics. Methods are being developed for exploring the mass spectrum 
and decay systematics of candidate superparticles.^ The key issue then will be whether 
their spins fit the pattern expected for SUSY. 

In this connection, a very interesting alternative hypothesis is that the new objects are 
indeed partners of the SM particles but with the same spins. This scenario is realized in 
the so-called universal extra dimension (UED) model, of the type proposed in refs. [4]. It is 
remarkable that UED and SUSY models could look very similar in collider experiments. In 
the former, the Kaluza-Klein (KK) excitations of the SM particles can carry a KK-parity 
analogous to the R-parity of the latter, implying pair production of the first KK-excitations 
and a stable lightest KK-particle (LKP), usually the first KK-excitation of the photon. If 
the energy of the LHC is insufficient to produce higher KK-excitations, then each SM 
particle would be found to have a single KK-partner with the same spin but a higher mass 
related to the size of the extra dimension. 

It is not our intention here to advocate or criticise UED models as a viable alternative 
to SUSY, but simply to discuss the extent to which they could be distinguished from 
SUSY models with the same mass spectra and decay systematics. Admittedly certain mass 
spectra would already be suggestive of one model or the other: for example, the spectrum 
of KK-excitations is quasi-degenerate if the extra dimension is simply compactified on an 
S1/Z2 orbifold [5]. The degeneracy is broken only by zero- mode (SM) masses, volume- 
suppressed boundary terms and loop corrections with a relatively low cutoff. In SUSY 
models one usually assumes that soft SUSY-breaking terms satisfy universal boundary 
conditions at a much higher scale, so that such degeneracy at the weak scale would be 
unnatural. However, the number of unknowns and arbitrary assumptions in each case is so 
great that the exclusion of cither class of model in favour of the other would only be truly 
convincing if their spin structures could be distinguished. 

In the present paper we assume that a particular chain of decays that is common in 
both SUSY and UED models, starting from a squark or a KK-excited quark respectively, 
has been identified and that all the masses of the new particles involved in it are known. 
We then study the extent to which decay correlations, manifest in the invariant mass 
distributions of combinations of observable decay products, would enable one to distinguish 
between the SUSY and UED spin assignments of the new particles. Our work is thus an 
extension of that described in ref. [1] (see also [6]), where the SUSY decay correlations 
were compared with uncorrelated phase space.^ 

In the following section we present the decay chains to be considered, and in section 3 
we recall the SUSY correlations investigated earlier. In section 4 we present our new results 
on the corresponding UED decay correlations. We derive simple analytical formulae for 
the correlation coefficients in terms of the masses in the decay chain, which should be of 
general use whatever the mass spectrum might turn out to be. We show graphical results 
for two specific mass scenarios, one considered more probable in SUSY and the other in 
UED. In both mass scenarios we compare the correlations predicted by the SUSY and UED 
spin assignments. 

'^See for example refs. [2, 3]. 

^Distinguishing between SUSY and UED at future e"'"e~ colliders has been studied in refs. [7, 8]. 
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As was emphasised by Barr [1], the obscrvabihty of interesting correlations depends 
cruciahy on the fact that the LHC is a proton-proton cohider, so that squarks/KK-quarks 
are produced somewhat more copiously than their antiparticles. To quantify this effect, 
we need to know the direct and indirect production cross sections of KK-quarks and KK- 
antiquarks. We have therefore computed the lowest-order two-parton to two-KK-parton 
matrix elements, which are expected to dominate the production of these particles. Our 
results, which differ somewhat from those presented in ref. [9],^ are discussed in section 5 
and listed in appendix B. 

Using our results on the UED production matrix elements and decay correlations, 
together with the decay branching ratios suggested in ref. [5], we have included a full 
simulation of the relevant UED processes in the HERWIG Monte Carlo event generator [10, 
11]. Since the corresponding SUSY processes, with full spin correlations, are already a well- 
established feature of HERWIG [12, 13], we are able in section 6 to present first detector-level 
results on distinguishing UED and SUSY spin correlations at the LHC. Our results and 
conclusions are summarized in section 7. 

2. Decay chains in SUSY and UED 




(a) (b) 

Figure 1: (a) SUSY and (b) UED decay chains considered here. 

The SUSY decay chain that we shall consider, which is the same as that studied in 
ref. [1], is shown in figure 1, together with the corresponding UED process. In both cases 
the visible decay products are a quark jet and a pair of opposite-sign same-flavour (OSSF) 
leptons with the same chirality. We suppose that the new particle masses have been 
measured, either by an edge analysis along the lines of refs. [2, 3] or some other means, and 
it remains to decide whether the decay angular distributions agree better with the SUSY 
or UED spin assignments. 

^An erratum to ref. [9] is in preparation (C. Macesanu, private communication). 
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The angular distributions depend on whether or not the chirahty of the slepton/KK- 
lepton is the same at that of the decaying squark/KK-quark.*^ For definiteness, we assume 
that the latter is left-handed, which is preferred in both of the models under consideration. 
We can then characterise the process by the chirality and charge of the "near" lepton, 
defined as shown in figure 1. Of course, we cannot distinguish experimentally between the 
"near" and "far" leptons, and so their contributions to any distribution will eventually 
have to be combined. However, in principle (in the zero-width approximation that we use) 
the processes with opposite "near" and "far" charge assignments are distinct. There are 
then two fundamental processes with different decay correlations, which (as in ref. [1]) we 
label 1 and 2: 

• Process 1: {q,l''<'^' = {qLJUt} oi' UlJI^^l} UlJrJr} or {qL,l^,l%}; 

• Process 2: {q, Z^-, Z^-} = {q^, it, l^} or {qL, l^, 1^ or {qL, l^, 1^} or {qL, 1^, l^}- 

3. Spin correlations in SUSY 

We first recapitulate from ref. [1] the angular distributions that are expected in the SUSY 
decay chain 1(a). The X2 spin one-half and its hclicity is the same as that of the 
quark, since the squark is a scalar. Therefore a near lepton with the same helicity as the 
quark (process 1) will be emitted preferentially at large values of the angle 9* between its 
direction and that of the quark in the X2 ^^st frame, with angular distribution (neglecting 
all SM particle masses) 

A pSUSY 1 

(l_cosr). (3.1) 

d cos 0* 2 

A near lepton with helicity opposite to the quark (process 2), on the other hand, will have 
angular distribution 

— 2-— = -(1 + cosr) . (3.2) 
dcosr 2^ ^ ^ ^ 

In terms of the ql^^^"^ invariant mass, 

{m^r? = 2|pHIp.I(1 - cosr) = l{m-nl^{l - cos9*) , (3.3) 
defining the rescaled invariant mass variable to be 

m = m^r/{m^r)^^ = sm{e*/2) (3.4) 

J pSUSY 

= 4m3 (3.5) 

am 

4m(l - m^) . (3.6) 



we therefore have 



and 



dm 



''We should emphasise that we use the term 'chirality' loosely here, since neither the spaxticles nor the 
KK-excitations concerned have definite handedness: what we mean is that they couple to SM particles of 
that chirality. 
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The slepton produced in the decay of the is a scalar particle, and so its decay is 
isotropic in its rest frame, and the near and far lepton directions are uncorrelated in that 
frame. There is nevertheless a weak correlation between the directions of the far lepton 
and the quark, due to the relative velocity of the slepton and the X2, which we discuss in 
more detail in appendix A. 

4. Spin correlations in UED 

4.1 Correlations in q* and Z* decays 

In the UED decay chain, the primary object has spin one-half and it decays to a quark 
and a vector boson, q* qZ*. The vector boson is either longitudinally or transversely 
polarised, with relative probabilities nig* and 2m^, , respectively. By angular momentum 
conservation, when the polarisation is transverse the helicity of the Z*, and of the near 
lepton, must be the same as that of the quark. The distribution of the angle 9* between 
the directions of the near lepton and the quark in the Z* rest frame is therefore of the form 

dPUEI^ _ 1 / dPo , , dP^ 



_(_ 2x ^ (4 1) 

dcos6'* 1 -F2a; Vdcos6'* dcos6'V ^'^ 

where x = m|,/m^* and P\ on the right-hand side represents the distribution for Z* 
helicity A. 

In the rest frame of the Z* , the decay of the longitudinally polarised state is forward- 
backward symmetric, with angular distribution 



(sin^0* + ycos^r) (4.2) 



dcosr 2{2 + y) 
where y = mf^/m^*. The transverse decay distribution for ibve helicity is 

-[{l±cos9*f + ysm'^e*] . (4.3) 



•^■^i ^ r/1 1 ___/)*^2 I 2/1*1 



dcosei* 4(2 -I- y)' 

The angular distributions for the two fundamental processes are therefore given by 

(jpUED o 

-T-^^ = 77 r [1 + X + T 2a;cos0* - (1 - x)(\ - y) cos^ 9*] . (4.4) 

dcose** 2{1 + 2x){2 + y) ^ \ )\ \ \ ) 

4.2 Quark + near lepton mass distribution 

Using eq. (4.4), the distribution of the rescaled qV^^^'^ invariant mass (3.4) is given in process 
1 by 

J pUED aff^ 

= (l + 2.)(2 + rt + - « + -y)^' - 4(1 - x)(l - rffti (4.5) 
and in process 2 by 



dm (l + 2x)(2 + y) 



[4x + y + 4{1 - 2x - y + xy)m'^ - 4{1 - x){l - y)m'^] . (4.6) 
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Notice that the fh distributions for the two UED processes become identical as x — > 0, 
i.e. for mz* ^ ruq*, since in that hmit the Z* polarisation is purely longitudinal. On the 
other hand when x ^ 1 and 7/ — > (m/* <C mz* ~ rUg*) the UED distributions become 
indistinguishable from those of the corresponding SUSY processes, eqs. (3.5,3.6). 

The features of the ql^'^^^ mass distributions can be illustrated by considering their 
forms for typical UED and SUSY patterns of new particle masses. A possible UED mass 
spectrum (from [5]) is given in table 1 for inverse radius = SOOGeV, cut-off A such 
that KR = 20 and rrih = 120 GeV. This model also assumes vanishing boundary terms 
at the cut-off scale A, and a vanishing boundary mass term for the Higgs mass, The 
lightest four left KK-quarks are degenerate in mass and are labelled here collectively as 
q^. Similarly the right and left KK-electrons and KK-muons are degenerate in mass and 
are labelled here as l*^ and l*^ respectively. This spectrum illustrates the feature of UED 
that the new particles have masses which are much closer to each other (usually within 
100 — 200 GeV) than in a SUSY spectrum based on high-scale universality. 

In the UED model we have eqs. (4.5) and (4.6) 
with X = m'^t./mg* = 0.803; the Z* decays preferen- 
tially to a left-handed excited lepton and so we use 
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Z* 




7* 

'■R 


7* 


501 


536 


598 


505 


515 
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m, 



: /m 



dPU^D 

dm 
dpu™ 

dm 



0.923, which yields 



0.727m + 2.577m'^ - 0.047m^ 



3.257m - 2.483m^ - 0.047m° 



Table 1: UED masses in GeV, for 
R-^ = 500GeV, AR = 20, rrih = 
120GeV, — and vanishing 

boundary terms at cut-ofF scale A. 



(4.7) 



These should be compared with the corresponding SUSY expressions (3.5) and (3.6), which 
are independent of the particle masses. 
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Figure 2: UED and SUSY distributions for (a) Process 1 and (b) Process 2 with respect to the 
rescaled gP®'^'' invariant mass, for the UED mass spectrum in table 1. Dotted: phase space. Dashed: 
SUSY. Solid/red: UED. 



The UED and SUSY angular distributions are plotted against each other for processes 
1 and 2 in figures 2(a) and 2(b) respectively. Since x = m^^/m^* is large in the typical 
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UED scenario, and the effect of y = m^t/m?z* weak at large x, the UED and SUSY 
distributions are similar. Therefore it will be difficult to verify the UED spin assignments 
if the spectrum is quasi-degenerate like that in table 1. 

The SUSY mass spectrum, on the other hand, does 
not naturally have the same near-degeneracy of neutrali- 
nos and squarks, and therefore the UED and SUSY an- 
gular distributions are more distinct. For illustration, we 
consider the MSSM Snowmass point SPS la [14], which 
has the mass spectrum shown in table 2. The decay 
X2 ^^R is preferred and therefore we use x = m^Q/m-^ 
0.653 for the comparative UED distributions, giving 



X? 




UL 






96 


177 


537 


143 


202 



Table 2: SUSY masses in GeV, for 
SPS point la. 



0.109 and y 



m~ /m~Q 
eRi X2 



dm 

dpUED 

dm 



1.213m + 3.108m'^ - 2.301m^ 



2.020m + 1.493m'^ - 2.301m^ 



(4.8) 
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Figure 3: UED and SUSY distributions for (a) Process 1 and (b) Process 2 with respect to the 
rescaled g/"''" invariant mass, for the SUSY mass spectrum in table 2. Dotted: phase space. 
Dashed: SUSY. SoUd/red: UED. 

The resulting mass distributions are compared in fig. 3. Owing to the small value 
of X, the UED predictions for the two processes are similar to each other, and different 
from the SUSY predictions. This gives some grounds for optimism that, if the spectrum 
is hierarchical, like that in table 2, then the SUSY spin assignments can be confirmed or 
ruled out in comparison with the UED assignments. 

4.3 Correlations in I* decay 

In the SUSY decay chain (figure la) , the slepton / is spinless and therefore it decays 
isotropically in its rest frame. In the UED case (figure lb), the spin of the KK lepton /* 
induces non-trivial correlations. Up to an overall constant, the full matrix elements for 
UED processes 1 and 2, as defined in section 2, take the form 

\M.\^ cx 2z(l 



z)Wi> + (1 - 2z)Wf 



(4.9) 
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where z = m^^/mf*, f represents the far lepton and, for / = T or /, 

Wi = (1 - X){2pz- ■ PnPZ* ■ Pi + rUz* Pn ■ Pi) - ^x{pn ■ PqPZ" ■ Pi + Pn ' PZ* Pq ' Pi) + 

_^ f . p^p^^ . fQj. process 1 

I Sx"^ Pn ■ Pq* Pq ■ pi for process 2 

where, as before, x = m\*/m^*, and n represents the near lepton. 

To specify the /* decay distribution, we define 9 as the angle between the near and far 
leptons in the I* rest frame, and (f) as the angle between the ql^^^ and dilepton planes, in 
the same frame. Then we find 

^3pUED o 

, ^^ , , , = — l + x + a;y^2xcos6'*-(l-x)(l-j/)cos2r 

dcos0*dcos0d(^ 47r(l + 2x)(2 + y) L ^ ^' 

1 — 2z ( 

-\^[l + x- xy^2xcos0* - (1 - x)(l + y) cos^ 6^*] cos 6' 



l + 2z 

-2[(1 -x)cosr ±a;]^sinrsin6lcos0}] . (4.11) 

4.4 Dilepton mass distribution 

The dilepton mass mu is simply related to the I* decay angle 6: 

ml = ^x\l - y)(l - z)il - cos9)ml. (4.12) 

and so, defining 

mi = 'mil/i'mii)raa.^ = sin(6'/2) (4.13) 

we have 

dP,f ° 4m,, 

The dilepton mass distribution is potentially a good indicator of UED spin correlations, 
because it is the same for processes 1 and 2 and relatively easy to measure. We see from 
eq. (4.14) that the deviation from the linear mass spectrum of phase space or SUSY is 
greatest when y and z are small, i.e., when m^* <C mi* <C mz*- On the other hand, the 
spin correlation vanishes when z = ^, that is, when m^* = mi*/V2, which is close to the 
situation at SPS point la. The typical UED scenario with quasi-degenerate masses also 
yields a small effect, with a mass spectrum proportional to mii{l — rhfi/5). Therefore a 
dilepton mass distribution differing significantly from phase space or SUSY would only be 
manifest in a UED model substantially different from those considered so far. 

4.5 Quark + far lepton mass distribution 

The quark + far lepton mass is a function of all the decay angles, 

1 



(mj')2 = l{l-x){l-z) (l + y)(l-cosrcos0) + 



+ (1 - y)(cos0* - cos6') - 2v^sin0*sin0cos(/) mi. (4.15) 
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Figure 4: UED and SUSY distributions with respect to the rescaled dilepton invariant mass, for 
(a) the UED and (b) the SUSY mass spectrum given above. Dotted: phase space. Dashed: SUSY. 
Sohd/red: UED. 



The maximum value occurs for 0* = and 9 = n, when 

(m'g7)max = - - z) ml* , 

and therefore the rescaled quark + far lepton mass is given by 



(4.16) 



1 



(l + y)(l-cos0*cos0) + 



+(l-2/)(cos6l*-cos6l)-2^sin6'*sin6'cos(/. . (4.17) 



The phase space for this quantity increases linearly up to the point rhqf = y^, then 
decreases logarithmically to zero. In the region fhqf < ^/y, the probability distribution has 
a polynomial form, as illustrated in figures 5 and 6 for the UED and SUSY mass scenarios 
respectively. Equations for the distribution are given in appendix A. 

We see that, as in the case of the dilepton mass distribution, the spin dependence of 
the quark + far lepton distribution is weak for both the mass spectra considered here. 

4.6 Observable quark-lepton correlations 

To proceed further, we must face the fact that the ql^'^^'^ and ql^^^ mass distributions are 
not experimentally observable. As pointed out in ref. [1], the best that can be done in 
reality is to measure the invariant masses of jet + lepton combinations. Assuming that the 
jet and lepton are indeed decay products from process 1 or 2, the jl^ mass distribution 
for a given lepton charge receives near-lepton contributions from the corresponding process 
and the charge conjugate of the other process, plus far-lepton contributions from the other 
process and the charge conjugate of the same process. Concentrating on the UED scenario 
of a left-handed KK-quark decaying to a left-handed KK-lepton, we have 



dP 



dm 



dP2 

^^near 



+ 



dPl 



dm 



far 



dPi 



+ 



dP2 



dm'^f^' ' dm^f 



(4.18) 
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Figure 5: UED and SUSY distributions for (a) Process 1 and (b) Process 2 with respect to the 
rescaled ql^'^'^ invariant mass, for the UED mass spectrum in table 1. Dotted: phase space. Dashed: 
SUSY. Solid/red: UED. 
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Figure 6: UED and SUSY distributions for (a) Process 1 and (b) Process 2 with respect to the 
rescsled ql^'" invariant mass, for the SUSY mass spectrum in table 2. Dotted: phase space. Dashed: 
SUSY. Solid/red: UED. 



where fq and /g are the quark and antiquark fractions in the selected event sample. Simi- 
larly 

dP ^ / dPi / dP2 dPi 

dmj7- l^dm^p^^ dm^f ) l^dm^f'^'' dm^f 

As will be discussed in section 5, for both the UED and SUSY scenarios we find fq ~ 0.7, 
fq ~ 0.3 at the LHC. The resulting observable jet+lepton mass distributions are then as 
depicted in figure 7, where again we have normalised to the maximum observable mass. 

Corresponding results for the SUSY mass scenario are shown in figure 8. Here the 
roles of processes 1 and 2 are interchanged in eqs. (4.18) and (4.19), since the decay now 
involves a right-handed slepton or KK-lepton. 
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Figure 7: UED and SUSY rescaled mass distributions for (a) jet + ^+ (b) jet + Z , for the UED 
mass spectrum in table 1. Dotted: phase space. Dashed: SUSY. Sohd/red: UED. 
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Figure 8: UED and SUSY rescaled mass distributions for (a) jet + (b) jet + / , for the SUSY 
mass spectrum in table 2. Dotted: phase space. Dashed: SUSY. Solid/red: UED. 



Figure 9 shows the resulting charge asymmetry 

^ ^ dP/dm^.;+ - dP/dm,^- 

dP/dmj;+ + dP/dm^,- ^ " ' 

We see that the UED and SUSY charge asymmetries are similar in form, but the latter is 
smaller by a factor of 2 - 4. 



5. Production cross sections 



As we saw in the previous section, the observability of a charge asymmetry sensitive to 
UED spin correlations depends critically on the difference between the production rates 
of KK-excited quarks and antiquarks, just as in SUSY it depends on the difference of 
squark and antisquark rates [1]. We have therefore computed the relevant lowest-order 
two-to-two scattering subprocess matrix elements and used them, together with the UED 



- 11 - 



A 



0.4 




, < , 1 , , , < 1 




0.4 






' ' ' ' / 
/- 
/ - 
/ - 
1 / ^ 


0.2 








0.2 


- 
















/ 

/ 


1 / ; 

\ / y 


0.0 








0.0 
























\\ 












-0.2 




\ 


\ \^ 


-0.2 














\ ^ 
\ 

\ 










-0.4 


^ , , , 1 , , , , 1 


. , . 1 . , . , 1 


, , , \ 


-0.4 









(a) m (b) m 

Figure 9: UED and SUSY charge asymmetries with respect to the jet + lepton rescaled invariant 
mass, for (a) the UED and (b) the SUSY mass spectrum given above. Dotted: phase space. Dashed: 
SUSY. Sohd/red: UED. 



branching ratios suggested in ref. [5], to estimate the UED production cross sections and 
the quantities /g and fq appearing in eqs. (4.18) and (4.19). 

Our expressions for the subprocess matrix elements are listed in appendix B. These 
results were obtained by including the Feynman rules for the effective four-dimensional 
theory in CompHEP [15]. They differ in several respects from those presented in ref. [9]. 
Details of the discrepancies are given in the appendix. Most importantly, we find a larger 
overall normalization. 

Our numerical results for the produc- 
tion cross sections at the LHC are pre- 
sented in table 3. These results were ob- 
tained from parton-level Monte Carlo sim- 
ulations of the production processes and 
decay chains, using the HERWIG event 
generator in SUSY mode with parton show- 
ering, hadronization and underlying event 
switched off. The HERWIG default (MRST 
leading-order [16]) parton distributions were 
used. For the UED simulations, the SUSY 
matrix element subroutine was replaced Table 3: Production cross sections (pb) in UED 
by a UED one and the SUSY particle data and SUSY models, with UED or SUSY masses, 
input file consisted of UED data based on 
ref. [5]. 

As a result of the more singular structure of the matrix elements and the extra helicity 
states, the UED production cross sections tend to be larger than those of the analogous 
SUSY processes for identical mass spectra, leading to an overall enhancement of the cross 
section, by a factor of about 8 for both the mass scenarios that we studied. Thus a SUSY- 
like signature (e.g. many jets and leptons plus large missing energy) with a cross section 



Masses 


Model 








/. 


UED 


UED 


249 


158 


83 


0.66 


UED 


SUSY 


28 


18 


9 


0.65 


SPS la 


UED 


480 


230 


102 


0.69 


SPS la 


SUSY 


55 


26 


11 


0.70 
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an order of magnitude larger than that predicted by SUSY models could be an initial 
indication of UED. 

The symbol ciaii represents the sum of cross section for all two-to-two subprocesses 
included for that model, while Uq* and Oq* represent the inclusive cross sections for direct 
and indirect production of the KK-quarks and -antiquarks (or squarks and antisquarks) 
initiating the decay chains in figure 1.^ Although the overall magnitude of the cross sections 
is different in the UED and SUSY models, we see that the KK-quark or squark fraction jq 
remains at about 70% for both models in both mass scenarios. 

6. Experimental observables 

To investigate the observability of the effects discussed above, we switched on the parton 
showering, hadronization and underlying event in our HERWIG simulations and applied 
cuts to approximate those that might be used to select new physics experimentally. For 
jet cuts we used the simple calorimeter simulation and cone jet finder GET JET [17], with 
cone size Ai? = 0.7. Our cuts were as follows: 

1. Missing transverse energy at least 100 GeV; 

2. At least four jets with transverse energies [Et] above 50 GeV; 

3. Sum of missing Et and four highest jet £^r's at least 400 GeV; 

4. Quark jet plus lepton invariant masses within the allowed range: mji± < {mqi)jnax- 

For cut 4 we cheated somewhat by selecting the jet that is nearest (in AR) to the true 
direction of the quark in the parton-level decay chain, and the leptons that do belong to 
that chain. Since we assume that all the new particle masses are already known, we expect 
that most ambiguities in the reconstruction of decay chains would in fact be resolved by 
invariant mass fits. 

Of course more sophisticated cuts and detector simulations could be applied to suit 
particular experiments, but this analysis should at least give an indication of whether 
further efforts in that direction are warranted. 

We included the UED spin correlations by generating decays according to phase space 
and reweighting events according to the invariant expression (4.9), evaluated at the parton 
level. To improve efficiency, we forced the decay Z* ^ 1*1 (l = e or n), and correspondingly 
X2 — > II, thereby enhancing the yield of the desired decay chains by factors of 3 and 8 in the 
UED and SUSY mass scenarios, respectively. To have comparable samples of each type, 
we generated 2 x 10^ chains of each type in figure 1 (including their charge conjugates) 
for each mass scenario. However, because of the different production cross sections and 
branching ratios in UED and SUSY, this corresponds to varying effective luminosities, as 
summarised in table 4. 

Also shown in table 4 are the fractions of events with the desired decay chains that 
survive the cuts listed above. For the SUSY mass scenario, the effect of these cuts is 

^Top-flavoured KK-quarks or squarks were excluded as their decays would have a different signature. 
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Masses 


Model 




Cut 1 


1+2 


1+2+3 


1+2+3+4 


UED 


UED 


7 


0.52 


0.11 


0.11 


0.05 




CTTCV 

dUd Y 


DO 


U.oo 


O.lo 


n 1 o 
U.iz 


O.Ob 


SPS la 


UED 


14 


0.86 


0.56 


0.56 


0.55 


SPS la 


SUSY 


131 


0.86 


0.54 


0.54 


0.53 



Table 4: Effective luminosities (fb ^) of our sample of 2 x 10^ decay chains, and fractions surviving 
the cuts listed above. 

not great. However, in the UED mass scenario the near-degeneracy of the mass spectrum 
means that the quark jet is likely to be relatively soft. It is therefore often misidentified or 
not found, with the result that few events survive the jet cuts. 

Figures 10 and 11 show the reconstructed jet plus lepton mass distributions for the 
UED and SUSY mass scenarios, respectively. As we saw in section 4, the expected quasi- 
degeneracy of the UED mass spectrum reduces the spin correlations and makes it hard 
to distinguish UED from SUSY in that case, even at the parton level. The experimental 
problems mentioned above, due to the relatively soft quark jet, further reduce the sensitivity 
and statistics. Consequently the charge asymmetry, shown in figure 12(a), is unlikely to 
be observable in the case of a UED-like mass spectrum. 




0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 



(a) m (b) m 

Figure 10: Detector-level rescaled mass distributions for (a) jet + l'^ (b) jet + for the UED 
mass spectrum in table 1. Dashed: SUSY. Solid/red: UED. 

The results for the SUSY (SPS la) mass spectrum, figures 11 and 12(b), are more 
encouraging. The spin correlations are larger in this case and, apart from some resolution 
smearing, their effects are not greatly diminished at the detector level. Correspondingly 
the charge asymmetry remains visible and similar to that predicted at the parton level, 
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0.2 0.4 0.6 O.B 1 0.2 0.4 0.6 0.8 1 

(a) m (b) m 

Figure 11: Detector-level rescaled mass distributions for (a) jet + Z+ (b) jet + , for the SUSY 
mass spectrum in table 2. Dashed: SUSY. Solid/red: UED. 

except at very high and low masses, where the asymmetry is the ratio of two vanishing 
quantities. 




0.2 0.4 0.6 O.B 1 0.2 0.4 0.6 O.B 1 



(a) m (b) m 

Figure 12: Detector-level charge asymmetries with respect to the jet -I- lepton rescaled invariant 
mass, for the (a) UED and (b) SUSY mass spectra given above. Dashed: SUSY. Sohd/red: UED. 



7. Conclusions 

We have presented results of a comparative study of spin correlations in models with 
supersymmetry and universal extra dimensions. Complete results were obtained for a 
decay chain that is likely to be important if either model is relevant at LHC energies. 
The analytical expressions for two-particle invariant mass distributions in section 4 can be 
used to test the models for any combination of masses and chirality of the new particles 
involved in the decay chain. We presented numerical and graphical results for two particular 
mass scenarios: one UED-like and one SUSY- like (SPS la). In the former case the near- 
degeneracy of the mass spectrum of new particles would make it difficult to verify the spin 
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content of the model in this way. In SUSY models such degeneracy would be less likely 
and the prospects for distinguishing between SUSY and UED spin assignments are better. 

As an adjunct to our study of spin effects we rederived the production cross sections 
for KK-partons in UED models and found some differences from results in the literature. 
Due to the more singular matrix elements and extra helicity states, the cross sections are 
significantly larger than those for SUSY particles with the same mass spectrum. This could 
also serve as a means of discriminating between SUSY and UED. 
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Appendices 

A. Quark + far lepton correlation 

The rescaled quark + far lepton invariant mass rhqf is given in terms of the decay angles 
by eq. (4.17). The phase space distribution has the form 

dpPS „^ Inj/ , ^ ^ 

— — = -2mqf- for < mqf < -s/y 

dmqf l-y 

= -4mg/i^^^ for < mg/ < 1 . (A.l) 

In the SUSY decay chain the only non-trivial dependence is on the angle 9*. We find 
the following form for the rescaled mass distribution: for process 1 



^pSUSY ^ 4^ 



Afhqf (1 - y) 



^(l-y-Flny) for < m^/ < 



^""^■^ (l-m^^ + 21nmg/) for < m^/ < 1 (A.2) 



and for process 2 

dpSUSY 4- 



(7 f 

+ forO<m<j/<^ 



dfhqf (1 - yf 
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In the case of UED, there is non-trivial dependence on ah the decay angles. We find 
the following form for the rescaled mass distribution: for process 1 



12m 



r{(l -y)[4x-y + 2(2 + 3y- 2x{5 + y))z 



dpU^D _ 

dfhgf " (l + 2x)(2+y)(l + 2z)(l-2/)■ 
-4m2^(2 - 3x){l - 2z)] - [y{l - 2(4 + y)z) + 4x{2z - y{l - 4z)) 
+4m2^(l + y-x{2 + y))(l - 2z)] In y} for < fhgf < ^ 



(l + 2x)(2 + y)(l + 2z)(l-y)2 

-4mlf{l - x)(l - z)] - 2[y{l - 2(4 + y)z) + Ax{2z - y{l - 4z)) 
+4m^^(l + y-x{2 + y)){l - 2z)] Inrhqf} for ^/y < rhgf < 1 

and for process 2 

^pUED ^ i2mgf 



{(1 - m^gf){Ax{l + 2y-bz- 6yz) - 5y + 2(2 + 9y)z 

(A.4) 



dfhgf (l + 2x)(2 + y)(l + 2z)(l-y)- 



- y)[-y + 2(2 + 2x(l - y) + 3y)z 



-4m2^(2 - x)(l - 2z)] - [y(l - 2(4 + y)z) 

+4m2 (1 + {l-x)y){l - 22;)]lny} for < m^/ < ^ 



12m, 



qf 



(l + 2x)(2 + y)(l + 2z)(l-y)2 



{(1 - mg/)[4(l + x)z - y{5 - 18z + 8xz) 



-4fhlf{l - x)(l - z)] - 2[y(l - 2(4 + y)z) 

+4m^^(l + (1 - x)y){l - 2z)] Infhgf} for < m^/ < 1 . 



(A.5) 



B. UED production cross sections 

We neglect all Standard Model particle masses and work at tree level, so that the n*'* 
excited KK state of each particle has mass n/R. In practice we ignore n > 1 excitations. 
We use here the notation of ref. [9]: ql or represent the first KK-cxcitations of the 
5D fields whose zero-modes are the left-handed doublet and right-handed singlet quarks 
respectively. Throughout, g| represents either q* or q^, Mi = 1/R, s,t,u are the usual 
Mandelstam variables, ts = t — Ml and ua = u — Ml. Only (B.l), (B.2) and (B.6) agree 
in form with ref. [9]. In addition we find an extra factor of 16 in the overall normalisation. 



2Ml tl + ul 



(B.l) 



at 



^Ml{- + ^-- 



23 2s^ 



+ 



8s 



6*3 8tl 



(B.2) 



9s_ 
27 



ul tl 



t3U4 
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+2 (3%+3^ + 4-^ 

V Ui ti t3U4 



(B.3) 



-4 



t^Ui \6t3n4 8 



s KfSt^Ui 8 



17 3t3M4 



Qt^Ui 24 4s2 



-9t 



+ 



3 [12t^ t|u4 



+ 



llsu4 5ti: 



6ti 



+ 



4Mi^ +5 + 4-2 +8- 



9 



-Mf ^ + 



^3J 



(B.4) 

(B.5) 
(B.6) 
(B.7) 



Mt 



2s' 



4s \ „ s 



+ 2 



+2 2 
13 U4 



t3n4 



+ 5 



(B.8) 



is 



(B.9) 



which is the same result as for qq — > qlql,qlql and qq' — > gi^i* as these are given by the 
same i-channel diagram. 



nil Ua\ 1 2ui 



which is the same result as for qq' qlq'*-^ 



(B.IO) 



9*9nf 



+1 + 



4 5 "h ^3 "h ^4 



3M 



(s^ + + ttl 

4s2t2^i2 



2£^+ti+^x| 
st3'U4 
t3n4 



(B.ll) 



9*9lf 



27 



4s3 57s 

+ 



108 



1,30,4 

+- 93 + 

^3114 



t3n4 s 
108^3^4 



(B.12) 



®We are grateful to Martyn Gigg and Peter Richardson for pointing out that the result for these processes 
had mistakenly been taken to be the same as equation (B.9) in an earlier version. 



- 18 - 



References 

[1] A. J. Barr, Determining the spin of supersymmetric particles at the LHC using lepton 
charge asymmetry, Phys. Lett. B 596 (2004) 205 [arXiv:hep-ph/0405052]. 

[2] B. C. AUanach, C. G. Lester, M. A. Parker and B. R. Webber, Measuring sparticle masses 
in non-universal string inspired models at the LHC, JHEP 0009 (2000) 004 
[arXiv:hep-ph/0007009] . 

[3] C. G. Lester, Model independent sparticle mass measurements at ATLAS, 
CERN-THESIS-2004-003 

[4] T. Appclquist, H. C. Cheng and B. A. Dobrcscu, Bounds on universal extra dimensions, 
Phys. Rev. D 64 (2001) 035002 [arXiv:hep-ph/0012100]. 

[5] H. C. Cheng, K. T. Matchev and M. Schmaltz, Radiative corrections to Kaluza-Klein 

masses, Phys. Rev. D 66 (2002) 036005 [arXiv:hep-ph/0204342]; Bosonic super symmetry? 
Getting fooled at the LHC, ibid. 66 (() 2002) 056006 [arXiv:hep-ph/0205314]. 

[6] T. Goto, K. Kawagoe and M. M. Nojiri, Study of the slepton non-universality at the CERN 
Large Hadron Collider, Phys. Rev. D 70, 075016 (2004) [Erratum-ibid. D 71, 059902 (2005)] 

[arXiv:hcp-ph/0406317]. 

[7] G. Bhattacharyya, P. Dey, A. Kundu and A. Raychaudhuri, Probing universal extra 
dimension at the International Linear Collider, arXiv:hep-ph/0502031. 

[8] M. BattagUa, A. Datta, A. De Roeck, K. Kong and K. T. Matchev, Contrasting 

supersymmetry and universal extra dimensions at the CLIC multi-TeV e+ e- collider, 
arXiv:hep-ph/0502041 . 

[9] C. Macesanu, C. D. McMuUen and S. Nandi, Collider im,plications of universal extra 
dimensions, Phys. Rev. D 66 (2002) 015009 [arXiv:hep-ph/0201300]. 

[10] G. Corcella, I.G. Knowles, G. Marchesini, S. Moretti, K. Odagiri, P. Richardson, 

M. H. Seymour and B. R. Webber, HERWIG 6: An event generator for hadron emission 
reactions with interfering gluons (including supersymmetric processes), JHEP 0101 (2001) 
010 [arXiv:hep-ph/0011363]. 

[11] G. Corcella et al, HERWIG 6.5 release note, arXiv:hep-ph/0210213. 

[12] P. Richardson, Spin correlations in Monte Carlo simulations, JHEP 0111 (2001) 029 

[arXiv:hep-ph/0110108]. 

[13] S. Moretti, K. Odagiri, P. Richardson, M. H. Seymour and B. R. Webber, Implementation 
of supersymmetric processes in the HERWIG event generator, JHEP 0204, 028 (2002) 
[arXiv:hep-ph/0204123]. 

[14] B. C. AUanach et al.. The Snowmass points and slopes: Benchmarks for SUSY searches, in 
Proc. of the APS/DPF/DPB Summer Study on the Future of Particle Physics (Snowmass 
2001) ed. N. Graf, Eur. Phys. J. C 25 (2002) 113 [eConf C010630 (2001) P125] 
[arXiv:hep-ph/0202233] . 

[15] A. Pukhov et al., CompHEP: A package for evaluation of Feynman diagrams and integration 
over multi-particle phase space. User's manual for version 33, arXiv:hep-ph/9908288. 

[16] A. D. Martin, R. G. Roberts, W. J. Stirhng and R. S. Thorne, Schem,e dependence, leading 
order and higher twist studies of MRST partons, Phys. Lett. B 443 (1998) 301 
[arXiv:hep-ph/9808371]. 



-19- 



[17] F. Paige, private communication. 



-20- 



